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ABSTRACT 

We present the detection of weak gravitational lensing in the field of the radio galaxy 
3C324 [z = 1.206) using deep HST imaging. From an analysis of the shapes of faint 
R = 24.5-27.5 galaxies in the field we measure a weak, coherent distortion centered 
close to the radio source. This shear field most likely arises from gravitational lensing 
of distant field galaxies by a foreground mass concentration. In the light of previous 
observations of this region, which indicate the presence of a rich cluster around the radio 
source, we suggest that the most likely candidate for the lens is the cluster associated 
with the radio galaxy at z = 1.2. If so, this is the most distant cluster to have been 
detected by weak shear observations. Such a statement has two important consequences. 
Firstly, it shows that massive, collapsed structures exist in the high redshift Universe, 
and secondly that a significant fraction of the R = 24.5-27.5 field galaxy population 
lies beyond z ~ 1.2. 

Subject headings: cosmology: observations - galaxies: evolution - clusters: general - 
gravitational lensing. 
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1. Introduction 

Observations of weak gravitational lensing are a 
relatively new tool to study the distribution of mass 
in the distant Universe. A growing number of papers 
employ this technique to investigate mass distribu- 
tions on a range of scales, from individual galaxies 
(Brainerd, Blandford fe Small 1995), through groups 
and poor clusters (Fort et al. 1995), rich clusters 
(Tyson, Valdes k Wenk 1990, Small et al. 1995a, 
Squires et al. 1995, Kneib et al. 1995) and onto larger 
scale structures (Mould et al. 1994). Lensing is partic- 
ularly useful in this regard because it does not require 
the mass in the lens to be virialised on the scales of 
interest. Lensing can also be used to constrain the 
redshifts of the very faint objects viewed through rich 
clusters by inverting simple geometrical models of the 
lens (Kneib et al. 1995). The primary limitation of 
this technique when applied to z < 0.5 clusters is the 
difficulty in obtaining accurate redshifts for the most 
distant objects, z > 1.5. This problem can be cir- 
cumvented if we can find a lens which is itself distant, 
z ~ 1 (c.f. Small, Ellis k Fitchett 1994). 

This letter presents one application of the weak 
lensing method to investigate the shear field in the 
vicinity of the distant radio galaxy 3C324 (z = 1.206) 
imaged with the Hubble Space Telescope (HST). HST 
is a unique instrument for this type of study provid- 
ing accurate shapes for faint galaxies, free of the sys- 
tematic effects which hamper ground-based applica- 
tions. We can thus expect to be able to measure weak 
shear signals at levels (< 1-3%) which are difficult to 
achieve on ground-based telescopes. 

We report a statistically significant shear pattern 
around 3C324. Previous observations (reviewed be- 
low) have indicated that 3C324 lies in a rich cluster of 
galaxies (Spinrad k Djorgovski 1984, Dickinson 1994, 
1995a). We suggest on the basis of these data that 
the likely source of the observed shear is the cluster 
at z = 1.2 associated with the radio source. We then 
present joint constraints on the properties of the lens- 
ing cluster and the distant field population which can 
be drawn from these observations. 

2. Observations, Analysis and Results 

2.1. Observations and Reduction 

The field around 3C324 was observed with WFPC2 
on-board HST for a total of 32 orbits in 5 visits be- 
tween 11 May and 15 June 1994. These deep images 



were obtained in order to study the morphologies of 
the faint, cluster galaxies (Dickinson et al. 1995). The 
final on-source integration time was 64.8 ksec, all in 
the F702W (R) filter. The individual images were 
acquired at the same roll angle with the same guide 
stars, but with small offsets between sets of observa- 
tions. Each WFPC2 exposure produces four separate 
frames, three for the Wide Field Camera (WFC) and 
one for the Planetary Camera (PC). These were mo- 
saiced together using a transformation which accounts 
for the geometric distortion of the WFPC2 optics. 
The mosaics were then aligned and combined with 
a clipping algorithm to reject cosmic ray events and 
other pixel defects. All reductions were carried out 
using standard IRAF/STSDAS routines. The stacked 
exposure is shown in Figure 1 (Plate XX). 

To catalog the faint galaxies within this frame we 
use the SExtractor image analysis package (Bertin 
1994). This software uses an isophotal area criterion 
to select objects detected above a fixed surface bright- 
ness threshold. It then provides photometry and 
shape information for all deblended objects. For this 
analysis we adopted a minimum area, after smooth- 
ing with a 0.3 arcsec diameter top-hat filter, of at 
least 0.12 arcsec'^ above the /jr = 26.4 mag arcsec"'^ 
isophote, equivalent to l.Ocr per pixel. To convert 
from instrumental magnitudes to R we take the con- 
version of Holtzmanet al. (1995) and use (V—R) ~ 0.5 
(Small et al. 1995b). 

The 3C324 field has a 5(7 detection limit oi R = 
27.7 within our 1.0 arcsec diameter photometry aper- 
ture, and contains ~ 1040 objects to this limit in the 
4.8 arcmin'^ of the 3 WFC chips. We discard both the 
lower sensitivity PC chip and also the borders of the 
WFC chips from our analysis. We apply three cuts 
to this catalog before using it in the lensing analy- 
sis. The first two cuts are on magnitude. We dis- 
card all objects fainter than R = 27.5 to remove ob- 
jects with inadequate signal-to-noise for measuring re- 
liable shapes. We also apply a bright magnitude cut 
to reduce foreground contamination, retaining only 
galaxies with R > 24.5, approximately 1 magnitude 
fainter than the brightest galaxies in z ~ 1 clus- 
ters. Finally, in order to remove objects which are 
too round to accurately determine their orientations, 
we remove all objects with ellipticities e < 0.05, where 
e = {a — b)/{a + b) where a and b are the maj or and mi- 
nor axis lengths of an object. The final catalog used 
in the lensing analysis contains ~ 700 objects with a 
surface density of ~ 145 arcmin"'^. 
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2.2. Analysis and Results 

To search for the effects of lensing we evaluate the 
complex shear: g = ee'^"^, where (j) is the position an- 
gle of the major axis of the object, for every object 
in the catalog (Kneib et al. 1995). These shear es- 
timates are then binned into independent 20" x 20" 
cells and plotted as a vector field over the image in 
Figure 1. Although a visual inspection of the im- 
age shows no good candidates for strongly lensed fea- 
tures, it is possible with some imagination to pick out 
a weak circular pattern in the shear vectors around 
the radio galaxy. To show this better we use the 
method of Kochanek (1990) to determine a center for 
the whole shear field. Each vector defines a proba- 
bility distribution for the likely center, in the form of 
two cones extending perpendicular to the shear direc- 
tion with opening angles determined from a combi- 
nation of the error in the shear direction, estimated 
from boot-strap resampling, and the shear strength 
at that point. The maximum of the combined prob- 
ability distributions from all the vectors then gives 
the best center for the whole shear field. The like- 
lihood distribution for the center is shown by prob- 
ability contours plotted on Figure 1. As suspected, 
the shear field is centered close to the radio galaxy, 
with its peak lying ~ 100/i~^ kpc^ to the south-west. 
In the remainder of the analysis we have adopted the 
radio source as the center of the lens. 

Having determined that a coherent shear field ex- 
ists around the radio source, we next wish to measure 
its strength and radial profile. We first calculate the 
average tangential shear [gi = — ecos(20), where 
is the angle between an object's major axis and the 
vector joining it to the lens center) within an annu- 
lus from r = 50-250/i~^ kpc (10-60 arcsec) around 
the radio source and obtain < gi >= 0.031 ± 0.010. 
Figure 2 plots gi averaged in circular annuli around 
3C324. The mean tangential shear drops by roughly a 
factor of 5 between 50-400/i~^ kpc, similar to the de- 
cline expected from lensing by a spherical, isothermal 
mass distribution. 

To assess possible systematic effects from the tele- 
scope optics, detectors, or the reduction method, we 
have applied the same analysis to a similarly deep 
HST image. This archival data (J. Westphal, P.I.) 
consists of deep V and /-band WFPC-2 images of the 



We take go = 0.5 and h = _ffo/100 km/sec/Mpc, which 
gives 4.30ft~^ kpc/arcsec at z = 1.206 



deepest field within the "Groth Strip," a blank field 
survey which does not target any known cluster. The 
data was processed and object catalogs were gener- 
ated in a similar manner as for the 3C324 dataset. Al- 
though the total exposure time is somewhat shorter, 
the galaxy surface density in this field is comparable 
to that available in the 3C324 observation. Within 
an annulus of 10-60 arcsec around the equivalent po- 
sition of 3C324 on the PC, we detect no significant 
tangential distortion (< gi >= -0.002 ± 0.012). We 
tentatively conclude that the observed shear in the 
3C324 field is neither instrumental in origin, nor a 
result of our reduction methods. 
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Figure 2. The radial dependency of the tangential 
shear strength (gi ) around 3C324. The error bars are la 
limits from boot-strap resampling of the data. We also 
plot a Qi (X. r" line with a = —1, showing the expected 
decline in signal if it arises from lensing by a spherical, 
isothermal mass distribution. This decline is formally con- 
sistent with the observations, which give a = —1.0 ± 0.2. 

3. Discussion 

We have demonstrated that a coherent shear pat- 
tern exists in the field of 3C324 and that this pattern 
is apparently not an artifact of the telescope optics. 
Moreover, the lensing pattern is centered close to the 
radio galaxy and shows a radial profile which is con- 
sistent with that seen from lensing by moderate and 
intermediate redshift clusters (Fort fe Mellier 1994). 
The available evidence therefore points towards the 
presence of a substantial mass concentration in the 
field of 3C324, and moreover this mass concentration 
lies in front of a sizeable fraction of the R = 24.5-27.5 
field population. We now address the question of the 
likely redshift of the lens. 

There is a marked excess of faint, red galaxies con- 
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centrated around 3C324 (Dickinson 1995a, 1995b). 
In the A'-band the overdensity relative to faint field 
number counts reaches a factor of ~6. Many of these 
galaxies have E/SO morphologies in the WFPC2 im- 
age and lie in a narrow color range with (R — K) « 6, 
too red for ordinary, unreddened galaxies at z < 1. 
Spectroscopic observations of these galaxies is cur- 
rently underway (Dickinson et al. 1995). In the re- 
mainder of our discussion we will assume that the 
shear field derives solely from a distant cluster, prob- 
ably associated with the radio source at z « 1.2. Here 
we rely upon the proximity of the radio galaxy to the 
center of the shear field to identify the 3C324 cluster 
as the lens. This is an important assumption: al- 
though the measured shear is statistically significant, 
it is also sufficiently weak that it could arise from a 
modest mass concentration at a lower redshift (c.f. 
Le Fevre et al. 1987, Fort et al. 1995). 

The observed shear strength depends jointly upon 
the distances of the background lensed galaxies and 
the mass within the lens. We will therefore investi- 
gate the allowed ranges for these two parameters by 
assuming plausible values for one parameter and de- 
termining the most likely value of the other. We start 
by examining the mass predictions. 

In view of the low signal strength and restricted 
field of view, complex modeling is not warranted. For 
simplicity, we will describe the lens by a spherical, 
singular isothermal mass profile with a single param- 
eter: the velocity dispersion, (Tci- We also take the 
following functional form (Efstathiou et al. 1990) to 
parametrise the redshift distribution, N[z), of the 
R = 24.5-27.5 field population: 

N{z) (X z2e-(^/^°)^ 

For our fiducial model, we adopt /3 = 5, for which 

< z >= 0.78zo. We focus on this model because it 
is a reasonable description of the ~ 25 field N{z) 
recently derived by Kneib et al. (1995). Those au- 
thors have inverted a well-constrained mass model 
of the rich cluster A2218 to estimate redshifts for 
~ 60 distorted arcs and arclets in the core. Their 
sample reached R ~ 25, at which point they obtain 

< z >~ 0.8, with ~ 20% of the population at z > 1. 
Extrapolating this result to the magnitude limit of the 
present study predicts <z>~ 1.2. Using our adopted 
form for N{z), Figure 3a plots the average tangential 
shear expected within our aperture as a function of 
(Tc;. The three curves show the dependence on the 
assumed mean redshift for the R = 24.5-27.5 field 



population, using < z >= 1.0, 1.2 and 1.5. These 
place 25%, 50% and 75% (respectively) of the faint 
galaxies at z > 1.2, and are so labeled. The high 
redshift of the lens, comparable to the mean redshift 
of the background sources, means that the mass es- 
timate is very sensitive to the exact choice for N{z). 
For the fiducial "no evolution" model with < z >~ 1.2, 
the observed shear corresponds to a velocity disper- 
sion ad ~ I3OOI230 ^™ sec~^, equivalent to a mass 
of M(r < 250/i-ikpc) ~ 3 x 10"/i-iM@. 




Figure 3 (a). The relationship between the average 
tangential shear, <gi>, measured in an annulus from 50- 
250h~^ kpc, and the lens mass, parameterised by ad- The 
functional form used for the galaxy redshift distribution is 
given in the text and we adopt /3 = 5 and an 80% efficiency 
for the shear measurement (Smail et al. in prep.). The 
three dashed curves show the predictions using different 
mean redshifts for the field population, labeled by the 
fraction of the population lying beyond z = 1.2. The solid 
line shows the measured shear, with ±1(7 bounds on this 
value shown by dotted lines. We derive ad = 1300^230 
sec~^ using the "no-evolution" model (50% with z > 1.2). 

It is clear from Figure 3a that for any plausible 
lens mass a substantial fraction (> 25%) of the faint 
field population must lie beyond z ~ 1.2. To fur- 
ther quantify this statement, we use an independent 
measure of the cluster mass and ask what is the shal- 
lowest faint galaxy N{z) which is consistent with the 
observed shear signal. Here, we make use of ROSAT 
observations of this field (Dickinson fe Mushotzky 
1995). An X-ray source coincident with 3C324 was 
initially detected in a 15.4 ksec ROSAT PSPC expo- 
sure, and was shown to be spatially extended using 
a deep 72.1 ksec integration with the HRI. Assum- 
ing that the entire flux arises from a hot intra-cluster 
plasma at the radio galaxy redshift with a temper- 
ature of 5 keV, metallicity 0.5 x solar, and cosmic 
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element abundance ratios, Dickinson fe Mushotzky 
(1995) determine a bolometric X-ray luminosity of 
Lx = (2.0 ± 0.4)/i-2 X 10^4 erg sec"!. 

For a nearby cluster with this luminosity, the local 
Lx-o'ci relationship (Edge fe Stewart 1990) would im- 
ply a rest-frame velocity dispersion of ad ~ 9OOI3QQ 
km sec~^. In view of possible differential evolution 
between the X-ray and optical properties of rich clus- 
ters observed to z ~ 0.5 (Castander et al. 1995), this 
may actually be a lower limit to the most likely value. 
Turning to Figure 3b, we find that a cluster with this 
(Tci would require a mean redshift < z >~ 1.51q'3 
for the R = 24.5-27.5 field population. Equivalently, 
at least 70 ± 20% of the population lies at z > 1.2, 
marginal greater than the 50% expected from the no- 
evolution model. Either way, it appears that at least 
half of the ~ 27 field population is at redshifts 
z > 1. 
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Figure 3 (b). The variation in <gi> with the mean 
redshift, < z >, of the faint field population. The three 
dashed curves trace the bounds on ad as estimated from 
the observed X-ray luminosity of the cluster and the lo- 
cal Lx-^d relationship: Ud ~ SOO^jqq km sec~^ . These 
translate into < z >= 1.5^q 3 for the galaxy population 
with 24.5 < _R < 27.5. 

This conclusion has important consequences for 
models of faint galaxy evolution. In particular, some 
recent models have proposed that the very faint field 
population is dominated by dwarf galaxies at z < 1. 
For example. Babul fe Ferguson (1995) consider a sce- 
nario where the faint counts are dominated by a mul- 
titude of "bursting dwarfs," whose formation is de- 
layed until z < 1. This model adequately reproduces 
observed galaxy counts and sizes, as well as their red- 
shift distributions to the limits of existing spectro- 
scopic surveys. At the fainter magnitudes probed by 



this lensing experiment, however, the model predicts 
much lower mean redshifts, with only ~10% of the 
24.5 < R < 27.5 galaxies lying at z > 1.2 (H. Fergu- 
son, private communication). 

4. Conclusions 

• We present a detection of weak shear in the field 
of the distant radio galaxy 3C324 at z = 1.206. We 
conclude that the shear arises from a mass concentra- 
tion in the direction of the radio galaxy, and appeal 
to other observations of this field to suggest that the 
shear is produced by a cluster associated with the ra- 
dio source at z = 1.2. 

• Simple modelling of the lensing cluster as an isother- 
mal mass distribution provides joint constraints on 
the redshift distribution of the faint field population 
to R = 27.5 and the amount of mass in the lens- 
ing cluster. Adopting the no-evolution N{z) for our 
magnitude limits we predict <z>~ 1.2. Taking this 
N{z) we estimate a velocity dispersion for the cluster 
of ad ~ I3OOI230 ^™ sec~^ from the observed shear 
signal, although this value is very sensitive to the ex- 
act choice of N{z). At a redshift of z ~ 1.2 collapsed 
structures of this mass are predicted to be extremely 
rare in standard cosmological models (Jing fe Fang 
1994). 

• Alternatively, we can fix the mass of the lens and de- 
termine the mean redshift of the faint galaxies in our 
sample. If the observed X-ray luminosity of 3C324 de- 
rives from a hot intra-cluster medium, and if the mass 
distribution is reasonably described by an isothermal 
sphere (an assumption supported by the radial depen- 
dence of the shear amplitude), then the local Lx^ad 
relationship would predict ad ~ 900 km sec~^. The 
observed shear strength would then require a mean 
redshift for the R = 24.5-27.5 field population of 
< z >~ 1.5lo 3, "with ~ 70% lying beyond z ~ 1.2. 
We therefore conclude that a large fraction (> 50%) 
of the ~ 27 field population must lie at z > 1. 
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Figures 

Plate 1. Combined HST WFPC2 image of 3C324 
taken through the F702W filter. The radio source 
lies on the PC chip and is marked. The total ex- 
posure time is 64.8 ksec and the effective resolution 
of this image is 0.1 arcsec. Overlayed on the image 
are the independent shear vectors determined from 
the shapes of the faint R = 24.5-27.5 galaxies. The 
dashed contours indicate the likelihood distribution 
for the best fitting center of the shear field. The scat- 
ter is estimated by repeatedly resampling the shear 
field and discarding 30% of the vectors. The contours 
are spaced by factors of 1.5 down from the peak. We 
show a vector indicating a 10% shear at the upper 
left. 
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